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Abstract 



Qh, We show that present experimental constraints on Bs ^ jj, jj, decay and the CDMS 

Oh! upper limit on the cold dark matter elastic scattering cross section already have significant 
^ ■ impact on the parameter space of the minimal supersymmetric extension of the Standard 
Model (MSSM) with non-universal supersymmetry-breaking scalar masses for the Higgs 
^ I multiplets (NUHM). The relaxation of scalar universality in the MSSM allows the possibility 
of a relatively light mass for the pseudoscalar Higgs boson. The present upper limit on 
Bg fi^ fi^ already excludes much of the scope for this possibility in the NUHM, in contrast 
to the constrained MSSM with universal scalar masses (CMSSM), where Bg decay 
does not exclude any ranges of parameters not already excluded by 6 ^ S7 decay. Cold dark 
matter scattering is also enhanced for small M^, but the impact of present upper limit on 
Bg — s> /x+ /i~ on the NUHM parameter space is in many cases greater than that of the CDMS 
scattering limit, particularly at large tan/5. 
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1 Introduction 



Many phenomenological analyses of the parameter space of the MSSM assume universahty 
for the soft supersymmetry-breaking scalar and gaugino masses, a theoretical framework 
often termed the constrained MSSM (CMSSM). However, this universality assumption is not 
necessarily supported by the effective supergravity models derived, for example, from string 
theory. On the other hand, the phenomenological suppression of flavour- changing neutral 
interactions suggests that squarks and sleptons with the same internal quantum numbers 
must be very nearly degenerate, at least for the supersymmetric partners of the first two 
generations, and there would be degeneracy (before renormalization) between squarks and 
sleptons in common GUT multiplets. However, there is no strong reason to suppose that 
the soft supersymmetry-breaking scalar masses mf,i = 1,2 of the Higgs multiplets should 
necessarily be the same as each other or the squarks and sleptons: mf = (1 + 6i)mQ with 
^1,2 7^ 0. These considerations motivate the phenomenological study of models with non- 
universal Higgs masses (NUHM) [1-6], as considered in this paper. 

The parameter space of the NUHM has two dimensions more than the CMSSM that are 
spanned by Si^2, allowing the Higgs supermultiplet mixing parameter n and the pseudoscalar 
Higgs mass uia to be treated as parameters that are free, apart from theoretical constraints 
such as vacuum stability up to the scale of grand unification: to this end, we impose the 
requirement that + /j,^ > at all renormalization scales below this GUT scale [4] . The 
phenomenological constraints on the NUHM provided by LEP constraints on the masses 
of the lightest supersymmetric Higgs boson ruh [7] and the lighter chargino have been 
considered, along with b ^ decay [8,9], the relic dark matter density flcDAth'^ [10] and 
(optionally) the anomalous magnetic moment of the muon, — 2 [11]. These constraints 
allow regions of NUHM parameter space in which uia is considerably smaller than its value 
in the GMSSM. 

The decay Bg fi^ is known to impose another interesting constraint on the param- 

eter spaces of models for physics beyond the Standard Model, such as the MSSM [12-15]. 
The Fermilab Tcvatron collider already has an interesting upper limit ~ 2 x 10~^ on the 
Bg decay branching ratio [16], and future runs of the Fermilab Tcvatron collider and 

the LHC are expected to increase significantly the experimental sensitivity to Bg — > yU"*" yU~ de- 
cay. However, a recent exploration of the present Bg fi~^ fi~ constraint in the CMSSM [17] 
found that its impact was limited by uncertainties in the theoretical relation of niA to the un- 
derlying CMSSM parameters, and provided no extensions of the (mi/2,mo) regions already 
excluded by 6 — > S7 decay, in particular. 
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We explore in this paper the current impact of this additional constraint on the NUHM, 
and consider also the potential impact of future improvements in the experimental sensitivity 
to Bg — > H'^ decay within the NUHM. The rate for Bs — > A*"*"//" may be enhanced in 
portions of the NUHM parameter space where Ma is smaller than in the CMSSM. We find 
that, consequently, significant regions of the NUHM parameter space at small tua and large 
tan (3 are already excluded by the present experimental upper limit on Bs — > decay. 
Likely improvements in sensitivity at the Fermilab Tevatron collider and the LHC will reach 
significant extra swathes of the NUHM parameter space. 

The elastic cold dark matter scattering cross section is also enhanced at small Ma, and 
another important constraint on the NUHM parameter space is placed by the upper limit 
on the spin-independent cold dark matter scattering cross section from the CDMS Collabo- 
ration [18]. This has only just begun to cut into the CMSSM parameter space [19], but does 
impact the NUHM parameter space, as also discussed in this paper. However, in many of 
the specific cases studied, the Bg — > jJi^ ji' constraint is stronger than the CDMS constraint. 

The structure of this paper is as follows. In Section 2 we recall briefly essential aspects 
of the theoretical calculation of Bs — > /x"^ decay and spin-independent dark matter scat- 
tering. Then, in Section 3 we explore various slices through the NUHM parameter space, 
displaying the impact of the present experimental upper limit on Bg — > //"^ decay. Fi- 
nally, in Section 4 we discuss the potential impact of future improvements in the experimental 
sensitivity to this decay. 

2 Review of the Calculations of Bs /x^ Decay 
and Spin-Independent Elastic % Scattering 

The branching ratio for the decay Bs — > n~ is given by 



X 
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where the one-loop corrected Wilson coefficients Cs^p are taken from [20] and Ca is defined 
in terms of Y{xt), following [21], as Ca = Y{xt)/ siv? 6w where 

y(a;,) = 0.99 f ^^V''' . (2) 
^ ^' Vl65GeV7 ^ ' 
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The function Y{xt) incorporates both leading [22] and next-to-leading order [23] QCD cor- 
rections, and mt{mt) is the running top-quark mass in the MS scheme. We assume here 
that the physical top quark mass is rrit — 172.7 ± 2.9 GeV [24]. 

The Wilson coefficients Cs,p receive four contributions in the context of MSSM, due to 
Higgs doublets [21], counter-terms, box and penguin diagram [25,26]^. We have considered 
all these one-loop corrections as well as the dominant NLO QCD corrections studied in [27] . 
In addition, we have included the flavour-changing gluino contribution [28,29]. The Wilson 
coefficients Cs,p have been multiplied by 1/(1 -|- e;,)^, where e;, incorporates the full flavour- 
independent supersymmetric one-loop corrections to the bottom-quark Yukawa coupling [30- 
32], that in principle are significant in the large-tan (5 regime [12,13]. Furthermore, it is known 
that the flavour- violating contributions arising from the Higgs and chargino couplings at the 
one-loop level result in effective one-loop corrected values for the Kobayashi-Maskawa (KM) 
matrix elements [26,33]. These corrections modify the Wilson coefficients involved in Eq. (1), 
as can be seen in Eqs. (6.35) and (6.36) in [34] or in Eq. (14) in [28]. We have included 
these flavour- violating effects as described in [28,34], taking into account the squark mixing 
effects. 

The counter-terms are mediated by A, H, h exchange, as seen in Eqs. (5.1) and (5.2) 
of [20], and dominate in the large-tan /3 limit, where the — > /J^'^ IJ'~ decay amplitude 
oc l/m\ and hence the decay rate oc l/m\. This underlines the potential of Bg — > /J^'^ IJ>~ 
decay for constraining models with large tan P and small m^. 

As already noted, in the NUHM uia may take values different from those required by 
the vacuum conditions in the CMSSM. Quite generally, the electroweak vacuum conditions 
in the MSSM may be written in the form 

2 _ ml -ml tan^ (3 + |m|(l - tan^ (3) + aJ/^ 
^ ^ tan2/3-l + A?) ' 

and 

m\{Q) = ml{Q) + ml{Q) + 2^,\Q) + A^(g), (4) 

and mi_2 = mi^2{fnz)-, where A^i, A^"^'^^ are loop corrections [32,35-38]. The exact forms of 
the radiative corrections to and are not needed for the discussion here, though we do 
note that the dominant contribution to A^''^^ at large tan /3 contains a term which is propor- 
tional to hi tan^ /3, whereas the dominant contribution to m\ contains terms proportional to 
hi tan 13 and h^ tan (3. The radiative corrections between the values of the quantities ji"^, ml 2 
at Q and the electroweak scale are well known. 

^See [20], where the full one- loop corrections have been calculated. 
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It is clear from (3) and (4) that departures of the input supersymmetry-breaking con- 
tributions to 2 from their universal values in the CMSSM, as permitted in the NUHM, 
induce corresponding changes in the allowed values of /i and m\, respectively. We evaluate 
all the relevant radiative corrections in our analysis: the sensitivity of to and nib was 
discussed extensively in [17], and we do not discuss the issue any further here. 

Spin-independent elastic x-nucleon scattering is controlled by the following effective four- 
fermion Lagrangian: 

^ = ^siXXQiQi, (5) 

which is to be summed over the quark flavours q, and the subscript i labels up-type quarks 
{i = 1) and down- type quarks (i = 2). The model-dependent coefficients a^i include terms 
oc l/m^^2 [39], where Hi 2 are the two scalar Higgs bosons in the MSSM, in which it is 
well known that the lighter one, H2, must have a mass ~ 120 GeV, whereas the heavier 
one. Hi, has a mass very similar to Ma- Hence the elastic cold dark matter scattering 
cross section also increases for small [5]. The magnitude of the cross section depends on 
hadronic matrix elements related to the 7r-nucleon E term, for which values between ~ 64 
and ~ 45 MeV are frequently quoted. In the estimates of the x-nuclcon scattering used here, 
we assume S = 64 MeV, which yields relatively large cross sections. Smaller regions of the 
NUHM parameter space would be excluded if we used a smaller estimate of S [19], so this 
assumption maximizes the possible impact of the CDMS constraint. 

3 Analysis of NUHM Parameter Planes 

In order to exemphfy the possible effects of the Bg — > and CDMS constraints, we dis- 

play some specific NUHM {/i, m^) planes for different values of tan mi/2 and mo, exhibiting 
the interplay of the different experimental, phenomenological and theoretical constraints. We 
first consider the case tan/3 = 10, mi/2 = 300 GeV and mo = 100 GeV, shown in panel (a) 
of Fig. 1. This value of tan/3 is towards the lower end of the range where we find generic 
GMSSM models satisfying all the constraints, for both signs of except for the g,j, — 2 
constraint that is satisfied only for positive /i. The most important constraints for > are 
those due to the chargino mass, shown as a vertical black dashed fine at low // ~ 100 GeV, 
and the GUT stabihty constraint at larger fj, ~ 650 GeV, shown as a near-vertical black 
dash-dotted fine that turns horizontal at low ~ 250 GeV. Regions between these fines 
are consistent with both constraints, but we also note a small excluded (brick-red) 'sugarloaf ' 
around // ~ 300 GeV that extends up to Ma ~ 300 GeV, where the LSP would have been the 
lighter stau, fi. For < 0, a large region is excluded by ft —> 57 decay, as shown by the green 
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Figure 1: Allowed regions in the {h^Ma) planes for mi/2 = 300 GeV and mo = 100 GeV, 
for (a) tan/3 = 10, (b) tan/3 = 40, (a) tan/3 = 50 and (d) tan/3 = 57. In each panel, 
the near-horizontal solid blue line is the contour where Ma = 2m^, and the turquoise strips 
are those where the relic neutralino LSP density falls within the range favoured by WMAP 
and other cosmological and astrophysical observations. The LEP chargino limit is shown 
as a dashed black line and the GUT stability constraint as a dot-dashed black line. The 
regions disallowed by b S'y are shaded green, and those disallowed because the fi or the 
Vf. would be the LSP are shaded brick-red and dark blue, respectively. The Higgs constraint 
is the largely horizonal red dot-dashed line. Gontours of the Bg /x^ n~ branching ratio 
are labelled correspondingly, with the current Tevatron limit the boldest black line, and the 
GDMS constraint is shown as a thick dashed grey line. In panel a), the — 2 constraint is 
satisfied when fi > 0. In the remaining three panels, the contribution to g^ — 2 is too large 
at the 2-a level. c 



(medium) shading here and in subsequent figures. The excluded regions where the z>e would 
be the LSP (or become tachyonic) are shaded (dark) blue: in this panel, they he beyond the 
GUT stabihty region. The red dash-dotted fine marks the LEP Higgs constraint ^, which 
in this case excludes the lower part of the region allowed by GUT stability, and requires 
Ma ^ 460 GeV. Finally, in this panel, regions with fj, > are allowed by 51^ — 2 ^, whereas 
regions with // < are disallowed at the 2-a level. 

Within the region allowed by the other constraints, we note that there is a near- vertical 
WMAP strip extending upwards. All of this WMAP strip is allowed comfortably by the 
current Bg — > //^ ii~ constraint, since the branching ratio for Bg — > /i" is very close to 
its Standard Model value ~ 3.9 x 10~^. It is dubious whether even the LHC would be able 
to distinguish the NUHM from the Standard Model for these values of tan/3, mi/2 and mo. 
Other regions with acceptable relic densities either have a small Higgs mass or violate the 
GUT constraint. These include the lower strip which extends upwards from Ma — 100 GeV 
and turns towards the horizontal at lower Ma, where rapid XX A,H annihilation becomes 
important when Ma ~ 2m^ (near-horizontal solid blue fine) and the WMAP strip in the 
upper left which is determined by neutralino-sneutrino co-annihilations. Also shown as a 
dashed grey line is the constraint imposed by the CDMS upper limit on spin-independent 
elastic cold dark matter scattering. Here and in the subsequent figures, in regions of the 
NUHM parameter space where the calculated rehc LSP density falls below the WMAP 
range, the cross section is rescaled by the factor ^x/^cdm, in order to compensate for the 
the fact that neutrahnos could provide only this fraction of the galactic halo. As a result, in 
this case, not only does the CDMS hmit not exclude any of the WMAP strip, it also does 
not exclude models surviving the LEP chargino constraint. 

When tan/3 = 40 ^, the stau constraint plays a much more important role, as seen in 
panel (b). It excludes all the parameter space above and to the left of a diagonal fine 
that meets the GUT stability constraint, which is essentially unchanged, at // ~ 700 GeV 
and Ma ~ 350 GeV. The z/g constraint appears at much larger and Ma- In this and 
subsequent panels of this figure, the LEP Higgs constraint excludes only a narrow strip 
below Ma ~ 120 GeV, whereas {g — 2)^ is incompatible with the measurement at the 2-a 
level. The only region allowed by the GUT stabihty constraint and by WMAP is then a 
short strip with fj, ~ 700 GeV, which lies above the Higgs constraint and below the sohd blue 
line where = 2m^. All of this strip is excluded by the current Bg — > /J''^ fJ'~ constraint 
^This is evaluated following the likelihood approach described in [40]. 

^We assume Sa^, where = {g^ — 2)/2, to be within the range 6.8 to 43.6 xlO~^° at the 2-a level. 
^For this and higher values of tan /3, we find consistent electroweak vacua mostly only for positive values 

of 
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(represented by bold solid line). We see also that the CDMS constraint excludes a region of 
the WMAP strip that is, however, already excluded by the stronger Bg — > /jl'^ constraint 
in this case. As already mentioned, the CDMS limit is rescaled when the calculated Qy^ 
falls below the WMAP range, which is responsible for the droop in the dashed grey line for 
fjL ~ 600 GeV. In this and subsequent figures, the CDMS hmit is also shown in regions where 
Qy^ exceeds the WMAP range, even though this region is disallowed by cosmology. 

In panel (c) for tan /3 = 50, the stau LSP constraint meets the GUT stabihty constraint 
at a somewhat higher values of Ma ~ 450 GeV. The region by WMAP and GUT stabihty is 
now bisected by the Ma — 2m^ line, with strips both above and below this line. The lower 
strip, where m^^ > Ma/2, is already excluded by the present Bg — > IJ>'^ IJ'~ constraint, and the 
upper strip lies only just beyond the present sensitivity. The lower WMAP strip would also 
be almost excluded by the CDMS constraint. Finally, in panel (d) for tan /3 = 57, which is 
close to the maximum value for which we find generic solutions to the electroweak vacuum 
conditions, we see that the stau LSP constraint now intersects the GUT stability constraint 
at Ma ~ 550 GeV, with the Og LSP constraint appearing only at larger values of There are 
again two strips allowed by WMAP, above and below the hne where Ma — 2m-x^. However, in 
this case, as well as the lower WMAP strip, a large portion of the upper strip is also already 
excluded by Bg — > /J''^ IJ'~, reflecting the greater power of this constraint as tan/3 increases. 
In this case, the CDMS constraint almost excludes the lower WMAP strip. 

Some qualitatively similar features are seen in Fig. 2, which displays analogous panels 
for the case mi/2 = 500 GeV, mo = 300 GeV. For tan /3 = 10, both signs of are equally 
possible, whereas we find no consistent electroweak vacuum in the polka-dotted region for 
tan P — 40, and no solutions with // < for the larger values of tan /3. In all the panels, 
the GUT stability constraint provides the right boundary of the allowed region at ~ 1100 
to 1350 GeV and, in panel (a) also a lower hmit on Ma ~ 350 GeV. For > 0, the LEP 
constraint on the chargino mass again supplies the left boundary in panel (a) for tan P — 10^, 
and partially in panel (b) for tan /3 = 40. The rest of the left boundary for > in panel (b), 
and the entire left boundaries in panels (c, d) for tan /3 = 50, 57, respectively, are provided 
by the stau LSP constraint. The bottom boundaries of the allowed regions in panels (b, c, 
d) are provided by the LEP Higgs constraint. 

As in panels (c, d) of Fig. 1, each panel features a pair of WMAP strips, above and 
below the Ma — 2m^ hne. When tan/3 = 10, the entire WMAP strips are allowed by 

''In panel (a), regions outside the near- vertical (pink) band at /x ~ 200 GeV are excluded by — 2. In 
the remaining three panels, all regions shown with /x > are allowed by — 2. In both panels (a) and (b) 
there are regions where /x < that are allowed by the other constraints, but disallowed by — 2. 



7 




Figure 2: As in Fig. 1, for the same values of tan (3 but for the choices mi/2 = 500 GeV 
and mo = 300 GeV. There is no electroweak symmetry breaking in the polka-dotted region of 
panel (b). In panel a), the vertical pink shaded region shows the region of g^ — 2 allowed at 
the 2-a level. In the remaining three panels, g^ — 2 is allowed when > 0. 
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Bg — > H''^ jJ-" I and sensitivity to the Standard Model prediction would be required to challenge 
any part of them. However, in the remaining panels, increasing parts of the lower WMAP 
strips are excluded by Bg — > //~ as tan (5 increases. However, in each case sensitivity to 
Bg — > ijf^ jJT below 10~^ would be required to explore all of the upper WMAP strip. The 
CDMS constraint lies outside the GUT stabihty region for tan /3 = 10. On the other hand, 
it excludes a somewhat larger part of the lower WMAP strip than does Bg — > jjf^ for 
tan P — 40, whereas the Bg — > constraint is stronger for tan P — 50 and 57. 

Both the stau LSP and GUT stability constraints weakened between Figs. 1 and 2, the 
culprit being the increased value of mo, which makes the stau heavier, and also makes the 
vacuum more likely to be stable. Therefore, in order to display more clearly the interplay of 
Bg — > with the hne where Ma — 2m^ and the pairs of upper and lower WMAP strips, 

we now consider a series of cases with relatively large mo = 1000 GeV, for which the stau 
and GUT stabihty constraints are irrelevant for large ranges of and Ma- Panels (a, b, c, d) 
of Fig. 3 show the cases mi/2 = 300, 500, 1000 and 1500 GeV, respectively, all for tan /3 = 57. 
The regions allowed by {g — 2)^ at the 2-a level are shown exphcitly in panels (b) and (c): 
at this level, the entire fi > region is allowed in panel (a) and disallowed in panel (d). We 
notice that the lines where Ma — 2m^ move upwards as mi/2 increases, reflecting the fact 
that m^ oc mi/2, approximately. On the other hand, the line representing the current upper 
bound on — > is relatively insensitive to both mi/2 and /i. For this reason, whereas 

Bg — > /i'^ ii~ already excludes the lower WMAP strip for the choice mi/2 = 300 GeV in 
panel (a) [compare also its impacts in panels (d) of Figs. 1 and 2], — > /i'^ ii~ is currently 
sensitive to only progressively smaller fractions of the lower WMAP strip as mi/2 increases 
in panels (b, c, d). By comparison, the CDMS constraint excludes only part of the lower 
WMAP strip in panel (a), but also a portion of the upper WMAP strip at small /i. The 
pecuhar shape of the CDMS curve in panel a) is caused by our scaling to the rehc density, 
which is particularly important when m^ Ma/2. In panels (b) and (c), CDMS excludes 
slightly more of the lower WMAP strip than does Bg — > jJi^ jJT , but this advantage would 
be removed if one adopted the lower value E = 45 MeV. The CDMS constraint has no 
significant impact in panel (d). 

To summarize our findings on the current impact of Bg /x+ ^~ in the (/x. Ma) planes: 
its importance increases with tan /3, and for large values it may exclude substantial parts of 
the WMAP strip where Ma < 2m^. The impact of Bg — > n'^ does not vary rapidly with 
mo, but is relatively less important as mi/2 increases. 

We now turn to some examples of parameter {Ma, tan /3) planes, shown in Fig. 4, which 
represent orthogonal projections of the NUHM parameter space that are often favoured in 



9 



tan P = 57 , 1111/2 = 300 , 111^ = 1000 



tan p = 57 , = 500 , hIq = 1000 




2000 



1000 



^ (GeV) 

tan P = 57 , = 1000 , hHq = 1000 



|a (GeV) 



tan P = 57 , mj/2 = 1500 , iHq = 1000 




-200 (B 



1000 



2000 



3000 



[I (GeV) 



|i (GeV) 



Figure 3: As in Fig. 1, hut for the choices mi/2 = (a) 300 GeV, (h) 500 GeV, (c) 1000 GeV 
and (d) 1500 GeV, all for tan (3 = 57 andniQ = 1000 GeV. The regions allowed by the — 2 
constraint are shaded pink (light grey) in panels (h) and (c). In panel (a), the region with 
fj, > is allowed by g^ — 2, whereas that in panel (d) is disallowed at the 2-a level. 
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Figure 4: Allowed regions in the {Ma, tan P) planes for (a) fi = 1000 GeV, mi/2 = 600 GeV 
and mo = 800 GeV and (h) = 800 GeV, mi/2 = 250 GeV and mo = 1000 GeV. The 
constraints are displayed in the same way as in in Fig. 1. 

analyses of MSSM Higgs phenomenology [41]. For convenience of comparison, our examples 
are taken from [4]. We see in panel (a) for /i = 200 GeV, mi/2 = 250 GeV and mo = 1000 GeV 
that 6 ^ S7 excludes regions at both large and small values of (M^, tan (3). The strip allowed 
by WMAP is relatively broad, with a portion to the right of the blue line where Ma = 2m^, 
and another part at low tan/5, but the latter is excluded by the LEP Higgs constraint. Only 
the parts of the plane with tan > 19 are allowed by — 2. In this case, only a small region 
at small Ma and large tan (3 is excluded by Bg —>■ fi^ fi^ . On the other hand, the GDMS 
experiment excludes all parts of the WMAP strip that are allowed hy — 2 ^. 

Turning now to panel (b) of Fig. 4 for ^ = 1000 GeV, mi/2 = 600 GeV and mo = 
800 GeV, we see that b ^ S'j now excludes only a region at small values of {Ma, tan P), 
and that the requirement of vacuum stability up to the GUT scale also excludes a region at 
Ma ^ 500 GeV. The WMAP strip is now much narrower than in panel (a), with portions on 
either side of the line where Ma = 2m^. The LEP Higgs constraint is weaker, excluding only 
tan/5 < 5, whereas the g^ — 'i constraint is stronger, excluding the region with tan/? < 35. 
In this case, neither the Bs fi^ ^~ nor the dark matter scattering constraints have any 
impact on the allowed portions of the WMAP strips at larger tan (5. 

As already mentioned, (M^,tan/?) planes are often considered in discussions of MSSM 
^This latter conclusion would, however, no longer hold if one adopted E = 45 MeV. 
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Higgs phenomenology. These two examples show that only small parts of these planes may 
be allowed by the various theoretical, phenomenological and cosmological constraints. In the 
particular examples studied, Bg — > is not yet an important constraint, whereas the 

search for astrophysical dark matter may be. Generic (MA,tan/3) regions may be allowed 
if one exploits the flexibility of the NUHM to vary /i, mi/2 and mo independently at fixed 
Ma and tan/3. We plan to return to a more detailed discussion of (Ma, tan /3) planes in the 
future. 

4 Summary and Prospects 

We have shown in this paper that the current Bs A*"*"/^" constraint already imposes 
significant constraints on the NUHM, excluding generic regions with small Ma and large tan P 
that would have been allowed by the other theoretical, phenomenological and cosmological 
constraints. The direct search for the scattering of cold dark matter also excludes some 
regions of the NUHM parameter space, but the Bg — > constraint is stronger in most 

of the cases we have studied. 

Experiments at the Tevatron and then the LHC are expected to increase greatly the 
sensitivity to Bg — > /J''^ IJ'~- Our analysis shows that, in many NUHM cases, this improved 
sensitivity would have good prospects for detecting Bg — > A*"*"/^". For example, sensitivity 
to BR(Bg jjL^ ~ 10~^ would give access to essentially all the (//, Ma) planes for 
mi/2 = 300 GeV, mo = 100 GeV and tan/3 ^ 40, as seen in Fig. 1. On the other hand, in 
the case of larger mi/2 = 500 GeV and mo = 300 GeV, shown in Fig. 2, the allowed regions 
of the WMAP strips extend to larger values of Ma that would require greater sensitivity to 
BR{Bg — > The same effect is seen even for the largest studied value of tan/3 = 57 

for several different choices of larger values of m 1/2 and mo, as seen in Fig. 3. Sensitivity to 
BR{Bg — > ~ 10~^ would not be sufficient to explore any new region of the {Ma, tan /3) 

plane for the choice fj, — 200 GeV, mi/2 = 250 GeV and mo = 1000 GeV, but would explore 
all the allowed region for fi — 1000 GeV, mi/2 = 600 GeV and mo = 800 GeV, as seen in 
panels (a) and (b) of Fig. 4, respectively. 

We conclude that Bg — > /x~ is already an important constraint on the NUHM param- 
eter space, and that, within this framework, it would have excellent prospects for a future 
detection of indirect effects of supersymmetry. 
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